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ABSTRACT 

M band spectra of two late-type T dwarfs, 2MASS J09373487+2931409, and 
Gliese 570D, confirm evidence from photometry that photospheric CO is present 
at abundance levels far in excess of those predicted from chemical equilibrium. 
These new and unambiguous detections of CO, together with an earlier spectro- 
scopic detection of CO in Gliese 229B and existing M band photometry of a large 
selection of T dwarfs, suggest that vertical mixing in the photosphere drives the 
CO abundance out of chemical equilibrium and is a common, and likely universal 
feature of mid-to-late type T dwarfs. The M band spectra allow determinations 
of the time scale of vertical mixing in the atmosphere of each object, the first such 
measurements of this important parameter in late T dwarfs. A detailed analy- 
sis of the spectral energy distribution of 2MASS J09373487+2931409 results in 
the following values for metallicity, temperature, surface gravity, and luminosity: 

[M/H] 0.3, T cS = 925-975 K, \ogg = 5.20-5.47, logL/L Q = -5.308±0.027. 

The age is 3-10 Gyr and the mass is in the range 45-69 Mj up i t er- 
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individual (Gliese 570D, 2MASS J09373487+2931409) 
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Introduction 



In the 1500-4000 K photospheres of late- type stars and early- mid L-type brown dwarfs 
carbon is predominantly locked up in carbon monoxide (CO), due to the much higher binding 
energy of CO than that of any other carbon-bearing molecule. At temperatures of ~1500 
K, however, methane (CH4) begins to be stable against collisional dissociation. As a brown 
dwarf cools below this temperature, t he overwhelming abunda nce of hydrogen slowly drives 
the carbon that was in CO into CH4 (iFegley fc Lodderslll996l ). This chemical change leads 
to a striking transformation of the infrared spectrum of the brown dwarf and marks the 



division between the L and T classificatio ns (IKirkpatrick et al.l Il999l ; iBurgasser et al.l 12002 
Geballe et allbooi burgasser 6^^120061 ). 



Absorption bands of CH4 together with those of H2O dominate the 1-2.5 /im spectra of T 
dwarfs, making the A = 2.3-2.5 /zm first overtone ban d of CO, which is ea sily observed in K to 
L dwarfs, virtually undetectable by spectral type T4 (IGeballe et al.ll2002l ). even if a significant 
amount of CO is still present. In contrast, the fundamental vibration- rotation band of CO 
is centered near 4.7/im, wher e methane is not an important absorber. Over a decade ago 
Noll. Geballe fc Marleyi (119971 ) discovered evidence of CO at 4.7 /im at an abundance roughly 
three orders of magnitude larger than expected from chemical equilibriu m, in the bright T7p 



brown dwarf Gliese 229B. The discov ery was promptly confirmed by lOppenheimer et al. 



(119981 ) . iNoll. Geballe & Marlevl (119971 ) found that the CO abundance in Gl 229B matched 
that expected at a temperature of 1250 K, roughly 400 K higher than the photospheric 
temperature. They interpreted the detection as evidence for vertical transport of CO from 
hotter in ternal atmospheric la yers to the surface, which was later supported by more detailed 
analysis (jSaumon et al.ll2000l ). 



This finding was not unexpected, as mixing is predicted to occur in the outer, ra- 
diative zones of brown dwarf atmospher e s, and the reaction(s) movin g carbon from CO 
to CH 4 are slow (jFegley k Lodderslll994! . Il996l : lLodders fc Feglevll2002h . Recently, the M 
band fluxes of numerous late T dwarfs have been measured to be up to one magnitude 
fainter than predictions b ased on their brightnesses in other wavebands and the assumptio n 
of equilibrium chemistry (jSaumon et al.l 120031 ; iGolimowski et al.l 120041 ; iPatten et al.l 120061 ) . 
This strongly indicates that the CO phenomenon seen in Gl 229B is common to most or 
all late T dwarfs. Here we describe spectroscopic observations of two late-type T dwarfs, 
2MASS 09373487+2931409 (hereafter, 2MASS 0937) and Gliese 570D, whose physical prop- 
erties differ significantly from Gl 229B, in order to test this hypothesis spectroscopically. 

2MASS 0937 and Gl 229B hay e nearly the same spectral classifications, T6p and T7.5, 
respectively (IBurgasser et al.ll2006l ). However, the K band flux of 2MASS 0937 is signifi- 
cantly depressed relative to Gl 229b and other T6-T7 dwarfs; e.g., its J — K color is 0.8 mag 
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bluer (IBurgasser et al.ll2002t iKnapp et al.l 120041 ) . This is believed to be due to the effect of 
greatly enhanced collision-induced absorption by H 2 (implying that 2MASS 0937 is consid- 
erably more massive than o t her T6-T7 dwarfs) and/or that it has unusually low metallicit y 
(IBurgasser et al.lhood I2OO3I : lBurrow^e^alll2002l:E^naDP et alJliool IBurgasser e^alll2006bh . 
The L — M color of 2MASS 0937 jGolimowski et all l2004h is close to that of Gl 229B, sug- 
gesting that its photosphere does contain an enhanced abundance of CO, but a spectroscopic 
observation is needed to verify that the this is indeed the cause. 

200_oj|_ is part o f a qua druple system whose other members are 



Gl 570D ( IBurgasser et al. 



low mass main sequence stars. ISaumon et al.l (120061 ) determined an age for it of 3-5 Gyr and 
a mass of 38-47 Mj upiter by comparing evolutionary models of brown dwarfs and 1-2.5 /j,m 
synthetic spectra of Gl 570D to the known luminosity and the observed spectrum. Compar- 



(Patten et al. 


2006: 


Leggett et al. 


2007a) 



5-15 /im spectrum of Gl 5 70D provides strong evidence for mixing also being responsible 



for the depletion of NH3 (ISaumon et al.l 120061 ) . In late T dwarfs the nitrogen chemistry 



is insensitive to the time scale of mixing in the radiative zone, because it is quenched in 
the convection zone deep in the atmosphere. On the other hand the carbon chemistry, and 
thus the abundance of CO, is very sensitive to the mixing time scale. Since Gl 570D is 
significantly cooler than Gl 229B and 2MASS 0937, observations of the 4.7/xm band of CO 
provide an opportunity to determine the mixing time scale for a different set of atmospheric 
parameters. 

This paper is structured as follows. Section 2 discusses the properties of the fundamental 
CO band, §3 presents the M band spectroscopic observations and data reduction, and §4 
covers the modeling and analysis of Gl 570D and 2MASS 0937, respectively. The last section 
summarizes and compares the results for the two objects. 



2. The Spectral Appearance of the Fundamental Band of CO 

In the photospheres of cool stars and "warm" brown dwarfs the fundamental (Av = 1) 
vibration-rotation band of CO, as viewed at low or medium resolution, is a broad and more 
or less structureless absorption stretching across the entire M band. This is the consequence 
of (1) several vibational levels v being populated, (2) each Av = 1 (i.e., 1-0, 2-1, 3-2, ...) 
band being successively shifted by ~0.06 /im (~25 cm -1 ) to longer wavelength, and (3) 
the spacing of vibration-rotation lines in each Av = 1 band being slightly different. If the 
rotational levels of only the ground vibrational state are populated, however, the absorption 
spectrum of CO observed at low or medium spectral resolution will exhibit a bump (due to 
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minimum absorption) around the center of the v — 1 band at 4.665 fj,m that separates 

the P and R branches of the band. As excited vibrational levels become more populated at 
higher temperatures, this absorption gap gradually disappears. 

The absorption gap is predicted to become noticeable at photospheric temperatures 
below 1500 K; that is, at about the same temperature at which equilibrium CO abundances 
are dropping rapidly in the photosphere of a cooling brown dwarf. Thus, for a mid-late type 
(T < 1000 K) T dwarf photosphere in chemical equilibrium, CO absorption lines will be very 
weak and no gap will be evident. If, however, the abundance of CO is considerably enhanced, 
then a broad absorption with a gap at the band center should be relatively prominent and 
detectable at moderate spectral resolution. The effec t is illustrated in early non-equilibrium 



model spectra of cool dwarfs of ISaumon et al.l (120031 ) 



H 2 also is an important absorber in the M band, but its lines, which are irregularly 
spaced, tend to be spread across the M band (although generally increasing in strength 
toward longer wavelengths). The H2O lines may partially mask the CO absorption gap at 
band center, but if the gap is prominent they should not entirely conceal it. 



3. Observations and Data Reduction 

M band spectra of 2MASS 0937 and Gliese 570D were accumulated over several nights 
during 2004-2005, at the Frederick C. Gillett Gemini North telescope, using the facility near 



infrared imager/spectrograph, (NIRI iHodapp et al. I (120031 )). which contains grisms for low 



resolution spectroscopy. The instrument was configured with a 0.75" slit and a grism, which 
together provided a resolving power of 460 (corresponding to AA=0.010 jttm) in the M band. 

An observing log is given in Table 1. The total exposure times were 6.1 hr for 2MASS 0937 
and 6.5 hr for Gl 570D. Observations were made in the standard stare-nod mode so that 
spectra were obtained alternately at two locations along the array, separated by 3". Coad- 
dition of subtracted pairs of frames then result in a final sky-subtracted image containing 
both a positive and a negative spectrum of the brown dwarf. Positioning of the target in the 
slit was checked using J-band imaging and was adjusted if necessary every 45-60 minutes. 
Typical adjustments were 0.1"; thus each brown dwarf probably was centered well within the 
slit at all times during most integrations. A telluric standard (a mid F dwarf) was observed 
just before or just after the brown dwarf, and at an airmass that closely matched that of 
the brown dwarf. The seeing, as judged from the intensity profiles along the slit of bright 
calibration stars, ranged from 0.5" to 0.7". 



To our knowledge, at M ~12.3 (as determined here) Gl 570D is the faintest astronomical 
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object to have its M-band spectrum successfully observed from the ground. Data reduction 
of M band spectra of faint sources such as Gl 570D and 2MASS 0937 is difficult, because not 
only is the residual background in the sky-subtracted coadded spectral image considerably 
brighter than the brown dwarf, but also the gradient in the residual background along the 
slit is steep compared to the signal produced by the spectrum of the faint source, even after 
~1 hr of integration. For these data it was necessary to extract sky rows adjacent on either 
side of each of the positive and negative spectra, in order to remove this gradient, prior to 
combining the positive and negative spectra. 

Wavelength calibrations of the spectra were obtained from the host of strong telluric 
lines evident in the spectra of telluric standards obtained just before and/or after the spectra 
of the brown dwarfs. The accuracies of the calibrations were better than 0.001 /im (3a). Flux 
calibrations utilized the 2MASS K$ magnitudes of the telluric standards and the very small 
color corrections to derive their M magnitudes. It was assumed that the spectra of the 
standards and the brown dwarfs suffered identical slit losses due to seeing and guiding. 
That assumption might be suspect, both because the brighter stars can be more accurately 
centered in the slit of the spectrograph and because during the longer integration times on the 
brown dwarfs it is more likely that guiding "drifts" occur and result in higher signal losses. 
Nevertheless, the average flux density in the 4.55-4.80 /xm spectrum of 2MASS 0937 in Fig. 1, 
5.0 ±0.1 x 10 -16 W m -2 /im" 1 , agrees to within several perc e nt wi th the M' (4.55-4.79 /im) 



magnitude of 11.74 ± 0.10) reported by iGolimowski et al.l (12004] ) . No M' photometry is 
available for Gl 570D. From Fig. 1 we estimate M' = 12.3±0.3 mag, where the unce rtainty 



is largely that of calibration. This value is consistent with the iPatten et al 



(20061) IRAC 



[4.5] photometry and values of [4.5]-M' for late T dwarfs as given by lLeggett et al.l (l2007bl . 
see their Fig. 2); the uncertainty in the color is several tenths of a magnitude. 

In Figure 1 the coadded, ratioed, and flux-calibrated spectra are binned int o 0.01 jum in- 



tervals. The spectra are qualita tively similar to the published spectra of Gl 229B (INoll. Geballe fc Marley 
19971 ; lOppenheimer et al.lll998l ). but have higher signal-to-noise ratios. Both have maximum 
flux density in the 4.64-4.70 /im region, as expected if photospheric CO depresses the spec- 
trum everywhere except near the 1-0 band center. The wavelengths of the 1-0 CO lines 
and their relative strengths, assuming optically thin lines and LTE conditions for T=750 K, 
are shown in the upper panel of Fig. 1 to illustrate the effect of the band center. Lines of 
H 2 also are expected to contribute spectral features. The model spectra (see below) show 
that the strongest H2O absorption feature in this region is at 4.68 /im, not far from the CO 
band center, and indeed the spectrum of each brown dwarf shows a pronounced dip at that 
wavelength. Thus the overall shapes of the spectra constitute unambiguous detections of 
CO in both 2MASS 0937 and Gl 570D. The effect of the band center is considerably more 
pronounced in the spectrum of Gl 570D than in 2MASS 0937. In the former the flux density 
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at the center is ~50% above the surrounding, whereas in the latter the contrast is roughly 
half of that. 



4. Modelling and Analysis 



We analyze the M band spectra with the goals of establishing the presence of CO on 
a quantitative basis and determining the time scale of vertical mixing in the atmosphere. 
The parallaxes of both objects are known and red, near-infrared, and mid-infrared Spitzer 
Space Telescope spectra are available, together sampling about 70% of the total flux emitted. 
The distance and spectra effectively constrain the luminosity, L, effectiv e temperature, T e ff , 
gravity, g, and metallicity. The Spitzer Infrared Sp ectrograph (IRS. iHouck et al.l 12004 ) 
spectra of T dwarfs later than ~ Tl show NH 3 features flushing et alJboO^ whose strengths 
are influenced by non-equilibrium nitrogen chemistry and thus by vertical mixing. Those 
features are not sensitive to the time scale of atmospheric mixing, however. Once other 
parameters have been constrained, the strength of the CO 4.7/im band can be used to 



find the mixing time scale. As in our previous work (iGeballe et al.l 1200 ll ; ISaumon et al. 
20061 ). we parametrize the vertical mixing time scale r m i x by the eddy diffusion coefficient 



K zz = H 2 /T T nw i where H is the mix ing scale height, assumed to be equal to the local pressure 



scale height ([Griffith fc Yellelll999h 



Our analysis is based on the atmospheric models and mixing scheme descr ibed in 
Saumon et al.l (120061 . 120071 ) and the evolution sequences of ISaumon fc Marleyl (120081 ). Since 
both Gl 570D and 2MASS 0937 are late T dwarfs, we assume that their atmospheres are 
cloudless, which is amply s upported by their JH K photom etry and the detaile d analysis 
of the T7p dwarf Gl 229B JSaumon etaJ boooh. G l 570D JSaumon et al.l l2006h . and the 
T8 dwarf 2MASS J0415195— 093506 (ISaumon et al.ll20071 ). The evolutionary sequences were 
computed with consistent surface boundary conditions provided by a grid of cloudless atmo- 
spheric models with the appropriate metallicities. 



4.1. Gliese 570D 



Because it is a companion to a well-studied main sequence star, Gl 570D has the most 
precisely determined physical parameters of any T dwarf to date. The K4 V primary provides 
an accurate distance and metallicity and a tight age constraint on the system. We obtain 
th e basic parameters fr om the bolometric luminosity L derived using the method developed 
by lGeballe et al.l (120011 ). Briefly, we use the well-sampled spectral energy distribution (SED) 
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and a combination of synthetic spectra and evolutionary sequences to obtain L. Given L, the 
evolutionary sequences give the range of possible models T eS (g) that are consistent with the 
age. A detailed comparison with the observed spectrum can further limit the possible range 
of solutions. This is a robust method to determine L, T e g and g for T dwarfs, as the results 
depend only weakly on the metallicity assumed for the models, the details of the synthetic 
spectra, and the degree of sampling of the SED. The method involves no assumption other 
than some confidence in t he models, which have been validated by comparison with data 
JSaumon et alJbood . l2007h . 



For Gl 570D. ISaumon et al.l fl2006f ) found T cff = 800-820 K, log#(cm/s 2 ) = 5.09-5.23, 
and log L/L Q = —5.525 to —5.528. The above ranges of values are not due to random 
uncertainties, because the parameters are correlated; e.g., for T cS = 800 K, log g = 5.09 and 
logL/ZvQ = —5.525. Fits of the spectrum favor the high-T e g, high-g limit of this range. 
The metallicity of the system is |Fe/H|=0. 9±0.04 (an average of s everal determinations : 
Thoren & FeltzineJ I2OO0I : ISantos et al1l2005l : IValenti fc Fisher! boosh . ISaumon et all J2006h 
used solar metallicity models to obtain the above parameters and we do likewise in our 
analysis of the M band spectrum. 

The 5-15 fim Spitzer spectrum shows strong ammonia features that can be accurately 
modeled using an NH3 abundance reduced by a factor of ~10 compared to that in chemi- 
cal equilibri um, as expected for non-equilibrium chemistr y driven by vertical mixing in the 
atmosphere (ILodders &: Fegleyl 120021 ; ISaumon et al.ll2006l ). We therefore expect that the M 
band spectrum will show a non-equilibrium abundance of CO. We have perfor med the analy- 



sis usi ng a baseline model computed with the most likely values determined by lSaumon et al. 
( 120061 ). which correspond to their model C: T e fj = 820 K, log g = 5.23, and an age of 5 Gyr 



(solar metallicity). The only free parameter left to fit the M band spectrum is the eddy 
diffusion coefficient K 77 , which determines the surface abundance of CO. 



4-1.1. Fits of the M band spectrum 

There are two approaches to fitting the spectrum shown in Fig. 1 with a set of synthetic 
spectra. The first method fits only the shape of the spectrum, i.e. the synthetic spectra 
are scaled by a constant multiplicative factor that is adjusted to minimize x 2 with the 
data. In this approach, the uncertainty in the fitting parameter (K zz ) is determined by 
the random noise in each data point of the observed spectrum. This is estimated at a = 
5.0 x 10~ 17 Wm -2 /im" 1 (Fig. 1). The second approach uses evolutionary sequences to obtain 
the radius of the brown dwarf and computes absolute model fluxes that are directly compared 
to the observed fluxes. This method simultaneously fits the shape of the spectrum and its 
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overall flux level and is subject to the systematic uncertainty in the flux calibration of the 
data, which is estimated at ±30%. We compute K zz and its uncertainty using both methods. 
With perfectly calibrated data and perfect models, both methods would give the same value 
of K zz . We find that the flux calibration uncertainty dominates the noise when fitting the 
absolute fluxes. The final source of uncertainty in the determination of K zz arises from the 
possible range of model parameters for the brown dwarf. For Gl 570D, where T e g is narrowly 
constrained, this last contribution is negligible. 

To account for the influence of noise on the interpretation of the spectrum, we generated 
1000 simulated observed spectra by adding random Gaussian noise with the above dispersion 
to the observed flux in each data point. Each of the simulated observed spectra was fitted 
with a synthetic spectrum computed with the baseline model parameters and \ogK zz (cm 2 /s) 
ranging between 2 and 8.5. The high resolution (A/AA ~ 16000) synthetic flux densities were 
converted to the resolution of the spectrum in Fig. 1 by integration over 0.01 /im wide bins. 
The best value of log K zz is obtained by minimizing 



y 2 = 1 y 



X r afr hs -(3fr d (K z 



a 

1 



where N is the number of points in the M band spectrum, /° bs , f™ d , and a are the ob- 
served flux, the absolute model flux, and the random noise in pixel i, respectively. The 
flux calibration uncertainty of ±30% enters as a (a = 0.7,1, 1.3) and /3 is the model flux 
renormalization. For the absolute flux fitting procedure f3 was set to unity. When fitting 
the shape of the spectrum only, (3 was adjusted to minimize \ 2 f° r each value of K zz . The 
resulting distribution of 1000 values of K zz gives an average value and its dispersion. The 
procedure was repeated for all three values of a (with (3 — 1) to quantify the effect of the 
uncertainty of the flux calibration. 

The general behavior of the \ 2 f° r these different assumptions is shown in Fig. 2. For a 
given value of K zz , the best fit, indicated by the lowest x 2 value, occurs when the model fluxes 
are scaled to the data, i.e. when only the shape of the observed spectrum is fitted (dotted 
curve). This curve has a broad minimum and our statistical sampling finds the minimum 



at log K zz (cm 2 /s) = 6.21 ±0.73. Equilibrium chemistry (K zz = 0), which ISaumon et al. 



(120061 ) found is strongly excluded on the basis of the NH 3 depletion revealed by the Spitzer 
spectrum, is also strongly excluded by the shape of the M band spectrum. This is readily 
seen in Fig. 3 which compares the equilibrium model spectrum, which shows no detectable 
CO, with the best-fitting non-equilibrium model spectrum. Furthermore, by increasing the 
statistical sampling by one order of magnitude to 10,000 simulated observed spectra, we found 
only three simulated spectra giving logK zz < 4 and none below 3.8. This illustrates the 
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negligible likelihood that the model using chemical equilibrium is consistent with the data. 
If we fit the data by considering the absolute flux calibration, we find narrower distributions 
of x 2 (solid curves on Fig. 2 that depend sensitively on the flux calibration of the data. 
For the nominal calibration, we get log K zz = 3.98 ± 0.12. Decreasing the absolute flux 
level by the estimated uncertainty of 30% we find \ogK zz = 5.23 ± 0.13. The systematic 
uncertainty in logK zz due to flux calibration is ~8 times larger than that due to random 
noise in the spectrum. Consistent solutions between the two fitting methods can be obtained 
if the calibrated flux density is reduced by a factor of 0.65 to 0.5, corresponding to a -l.ler 
to -1.7(T correction, which is just plausible. 



4-1.2. Photometry 



Photometric measurements of Gl 5 70D in the MKO V (IGolimowski et al.l 120041 ) and 
the Spitzer IRAC [3.6] and [4.5] bands (IPatten et al.l 120061 ) provide additional relevant in- 
formation. We generated synthetic fluxes by integrating the model spectra over the band- 
passe^], and multiplying by [Rj D) 2 , where R(T e g, g) is the radius obtained fro m our cloudless 
[M/H]=0 evolution sequences and D is the distance (jPerryman et al.l 119971 ) . The observed 



magnitudes were converted to fluxes by computing zero magnitude fluxes in each bandpass 
with a spectrum of Vega|^| or from the flux calibration of the IRAC instrumental A least 
squares fit of the three photometric points gives \ogK zz (cm 2 /s) = 4.49 ± 0.24 for our nom- 
inal model (T e ff = 820 K, log g = 5.23 and [M/H]=0). This best-fitting synthetic spectrum, 
its corresponding equilibrium spectrum and the photometry are shown in Fig. 4 along with 
the M band spectrum. The M band spectrum in the figure is scaled by a factor of 0.83 
relative to the flux-calibrated spectrum in order to minimize the residuals with the model 
spectrum. This is a shift of about half of the calibration uncertainty of the spectrum. The 
figure clearly shows the strong and broad CO band, centered at 4.7/xm and complete with 
the bump at band center, that appears when vertical transport drives the carbon chemistry 
far out of equilibrium. The spectrum outside of the 4.4-5.0 /im region is barely affected, 
however. 



1 MKO system: Tokunaga et al. ( 2002h : IRAC: [http://ssc.spitzer.caltech.edu/irac/spectraL response.html 
2 www. stsci.edu/instruments/observatory/cdbs/calspec 
3 http:/ /ssc. spitzer. caltech.edu/irac/calib/ 
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4-1.3. Summary: Optimal model 

Figure 4 shows that a cloudless model with T e g = 820 K, log g = 5.23, [M/H]=0 and 
logK zz (c m 2 /s) =4.5 agrees w ith (1) the optical, near infrared and the 5-15 jum Spitzer IRS 
spectrum ( Saumon et al. 20061 ). (2) the MKO L' and Spitzer IRAC [4.5] photometry, and (3) 



the flux level of the M band spectrum within its uncertainties. The IRAC [3.6] flux is too low, 
however. In addition log K zz = 6.21 ±0.73 provides a better fit to the shape of the spectrum. 
This latter value also fits all of the data except for the photometry. Within the various 
sources of uncertainty, these two results for log K zz are barely compatible. For instance, the 
(non- Gaussian) statistical distribution that led to the determination \ogK zz = 6.21 ± 0.73 
has only 1.6% of the sample falling within the range 4.5 ± 0.24. 

We believe that fitting the shape of the M band spectrum only is the more reliable way 
to determine the value of K zz . This method depends primarily on the CO/H2O abundance 
ratio in the photosphere and much less on modeling details While the flux level in the 
Spitzer 4.5 //m band is indicative of CO absorption, it does not specifically point to CO 
as the absorber and the filter is not well matched to the CO band. In addition, there is 
no photometric data point outside of and just to longer wavelength of the CO band, where 
absorption by H2O, which also occurs within the CO band, should dominate. 

The photometric measurements are much more precise than the M band spectrum 
calibration, however. They measure the spectral energy distribution in the vicinity of the 
4-5 fim peak, which is strongly affected by the 2/3 band of CH4 centered at 3.3 /j,m extending 
to ~ 4 fim, and by the v 2 band of H 2 on the long wavelength side of the peak. The CH 4 line 
list is known to be incomplete at the temperatures o f brown dwarf atmosp heres. The modest 



success in fitting the 3-4 /mi spectra of T dwarfs ([Stephens et al.l 120081 ) suggests that the 
inconsistency between the two solutions could be resolved in favor of the higher value of K zz 
by invoking a higher CH4 opacity. Indeed the incompleteness of the CH4 line list implies that 
the CH 4 opacity in our model is a lower limit to the actual opacity. Increasing the opacity 
in the red wing of the v 3 band of CH 4 would nudge all three broad band fluxes in a direction 
that would increase the value of the best fitting K zz . The flux removed from the increased 
CH4 absorption in the near infrared will reemerge primarily in low opacity windows. It is 
very likely that the shape of the 4-5 /zm peak will be altered by flux redistribution caused 
by a change in the CH 4 opacity, which could significantly change the value of K zz that best 
fits the photometry. 

Thus, we adopt \ogK zz (cm 2 /s) = 6.2 ± 0.7 as the most likely value for G l 570D, based 



on the favored nominal model parameters (model C of ISaumon et al.l (12006)). Repeating 



the ab ove analysis with the less favored but still possible values of model B of 



Saumon et al. 



(|2006|) (T eff = 800 K, \ogg = 5.09), gives log^cm 2 /s = 6.0 ± 0.7 which is essentially 
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identical within the uncertainty. The corresponding time scale of vertical mixing in the 
radiative region of the atmosphere of Gl 570D is ~ 0.1 to 3 hours. The complete infrared 
synthetic spectrum of our optimal model for Gl 570D is not shown b ecause it is essentially 



indistinguishable from that shown in Fig. 2 of ISaumon et al.l (120061 ) which was computed 
with \ogK zz = 2. The only region of the spectrum that is affected by the larger value of 
K zz is between 3.8 and 5.0 /im, which is shown in Figs. 3 and 4. 

For the optimal model the composition of the atmosphere as a function of depth for the 
most important opacity sources is shown in Fig. 5. The chemistry of nitrogen is quenched in 
the convection zone where r m i x becomes shorter than tn 2 , the time scale for conversion of N2 
to NH3. This occurs at logT ~ 3.38. In the upper atmosphere, the NH3 abundance is nearly 
a factor of ten lower than the equilibrium value. Because the eddy diffusion coefficient K zz 
determines the mixing time scale in the radiative zone, the quenching level of the nitrogen 
chemistry is unaffected by the choice of K zz ; rather it is determined by the convective mixing 
time scale. On the other hand, the carbon chemistry is quenched in the radiative zone around 
logT = 3.17 and its abundance depends sensitively on K zz because the equilibrium mole 
fraction of CO has a steep dependence on depth in the atmosphere in low-T e g models. In the 
upper part of the atmosphere the CO mole fraction is logXco = —4.10 ± 0.18. Finally, the 
4.7/im band of CO and the 10-11 /mi feature of NH 3 probe the same level of the atmosphere, 
because the lines in both spectral regions are formed at logT ~ 2.75-2.9. 



4.2. 2MASS 0937 



4-2.1. Determination ofT e g, gravity and metallicity 



Due to its very blue J — K color, 2MASS 0937 has long been recognized as hav 



2002 



ing low metallicity and /or high gravity (IBurgasser et al.l 120021 . 120031 ; iBurrows et al, 
Knapp et all2004 ; Burgasser et al. 2006b ). Previous esti mates give T p r = 725 -1000 K ( Golimowski et al 
20041) and 70 0-85 K f|Vrba et al-l hooll). and \ogg ~ 5.5 fcnapp et al.l 120041 ). More recently 



Burgasser et al.l (j2006bl ). using a calibrated set of spectral indices, found that 2MASS 0937 
is metal poor with [M/H] between -0.4 and -0.1. For [M/H]= -0.2, they obtain T cS = 780- 
840 K and logg = 5.3-5.5. With the wealth of information provided by the wide spectral 
coverage now available, we can determine these parameters more reliably. The observed 
spectrum covers the intervals 0.63-1.01 /im, 1.03-1.345 /im, 1.40-2.53 /tm and 5.13-15.35 /im. 
The integrated observed flux is 2.97 x 10~ 12 ergs -1 cm -2 , which represents about 70% of the 
total flux. The flux calibration uncertainties are ±3% in the red, ±5% in the near infrared, 



and ±3.7% for the Spitzer IRS spectrum. The distance is 6.141 ±0.146 pc (IVrba et al.ll2004l ). 
The resulting uncertainty in the bolometric luminosity is ±6.4%, or 0.027 dex, where we have 
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combined the various flux calibration uncertainties rather than treating them as statistically 
independent variations. 

We first applied our suite of solar metallicity models and found, as expected, that the 
model overestimate the K band flux by about 60%, even at the highest possible gravity 
of log g = 5.43. A better match of the near-infrared spectrum, and the K band flux in 
particular, thus requires a sub-solar metallicity. Applying our gr id of model atmosph eres, 
spectra, and evolutio n for [M/H]=— 0.3 to the method described in lGeballe et al.l (120011 ) and 
Saumon et al.l (120061 ). we obtain the parameters shown in Fig. 6 and given in Table 2 for 
four different ages ranging from 1 to 10 Gyr. We will see below that models corresponding 
to ages of 1 Gy and less do not fit the data satisfactorily. We find th at log L/L m ~ — 5.30 , 
the same luminosity as the T7.5 dwarf 2MASS J12171110— 0311131 JSaumon et al.l l2007h . 
The red and near infrared spectra of these four models are nearly identical, as can be seen in 
Fig. 7 where the two extreme cases, models A and D, are plotted. The lower gravity model 
matches the data best, but the four models reproduce the J and K peak fluxes equally well. 
Note that the fluxes in the model spectra in Fig. 7 are absolute and have not been scaled to 
the observed spectrum. The model spectra were computed from parameters determined with 
the integrated spectral flux and the distance to 2MASS 0937 as the only input. The excellent 
agreement implies that the metallicity of 2MASS 0937 is close to [M/H]= —0.3, the value used 
for the rest of the analysis. For comparison, the values obtained with the solar metallicity 
models are shifted by 10-15 K lower in T e s, ~ 0.05 dex lower in log g and < 0.02 dex lower 
in logL for a given age. These small variations, which are close to the uncertainties in the 
present determinations, demonstrate that a more accurate determination of the metallicity 
of 2MASS 0937 will not significantly affect the parameters given in Table 2. 

The range of parameters between models A and D cannot be further constrained on the 
basis of the near infrared spectrum. Generally, the K band flux (e.g., J — K) is a gravity 
indicator for a fixed T e g. The models in Table 2 are nearly constant in L. Higher gravities 
also imply higher T e ff, with compensatin g effects on the near infrared spectra (Fig. 7). As 
we had found in our analysis of Gl 570D (jSaumon et al.ll2006l ). the mid-infrared Spitzer IRS 
spectrum is more sensitive to variations in logg in the suite of T e s(g) solutions than is the 
near-infrared spectrum. We have compared each of the four models in Table 2 with the mid- 
infrared data, taking into acco unt the spectrum o f random pixel noise and the systematic 
±3.7% calibration uncertainty (IPatten et al.ll2006l ). Because there is a strong suspicion that 
the chemistry in the atmosphere of 2MASS 0937 is driven out of equilibrium by vertical 
mixing, we also consider out-of-equilibrium synthetic spectra computed with the parameters 
of models A through D. The mid-infrared sp ectrum is insensitive to the choice of the vertical 
mixing parameter K zz (jSaumon et al.l 120061 ) and we arbitrarily choose a value of 10 4 cm 2 /s. 
To constrain the set of T e s(g) solutions, we thus consider eight model spectra (the four 
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models of Table 2, with both K zz = and 10 4 cm 2 /s) and the observed 5-15 /im spectrum 
with its nominal calibration, the calibration modified by ±3.7%, and finally by fitting only 
the shape of the spectrum, i.e. without regard of the absolute flux level. 

The goodness of fit of a particular model is measured by the x 2 > which includes the 
random uncertainty in the flux of each pixel. The ensemble of x 2 f° r these models reveals 
a clear picture. In every case, the best of the equilibrium models (A-D) is a much worse 
fit than the best non-equilibrium model, at more than an 8a level of significance. The IRS 
spectrum thus rules out equilibrium chemistry for nitrogen in 2MASS 0937. For the non- 
equilibrium case, models A through D fit the shape of the spectrum equally well. If absolute 
fluxes are considered, for all choices of calibration (within ±3.7%) the lowest gravity model 
(A) is ruled out at > 22a. Models B through D can provide equally good fits depending on 
the choice of absolute calibration. For the nominal calibration, the best fit is obtained with 
model C. We conclude that the parameters of 2MASS 0937 fall within the range between 
models B and D in Table 2 and that vertical transport occurs in its atmosphere, resulting in 
a depletion of NH3. The synthetic spectrum for non-equilibrium model C is shown in Fig. 8 
as well as the best equilibrium model (model B). The non-equilibrium model agrees very 
well with the data, although the fit is not perfect. 

To summarize, combining the information in the red, near infrared, and mid-infrared 
spectra with the distance to 2 MASS 0937 we find that T eff = 925 to 975 K, log g = 5.19 
to 5.47, log L/L & = -5.31 ± 0.027, the metallicity is close to [M/H]= -0.3, and the age 
is between ~3 and 10 Gyr. These values confirm that 2MASS 0937 is a metal-poor brown 
d warf with a fairly high gr avity. The metallici t y and gravity agree with the determinations 
of iBurgasser et al.l (j2006bl ) and iKnapp et al. d2004h but T e ff is su bstantially higher than 
previous determinations, except for that of iGolimowski et al.l ( 120041 ) . 



4-2.2. Fits of the M band spectrum 

We fit the M band spectrum of 2MASS 0937 to determine the value of the eddy diffusion 
coefficient K zz using the procedure described in §4.1.1. In this case, the random noise in 
each pixel of the observed spectrum is estimated at a = 4.2 x 10 _17 Wm _2 //m _1 and the 
systematic uncertainty in the flux calibration is estimated at ±30% although the agreement 
with previous M' photometry is much better than this. We consider the three models B 
though D in this analysis. By fitting the shape of the M band spectrum only, we find that 
\ogK zz (cm 2 /s) = 4.3 ± 0.3 for model C. The value increases slowly with logg among the 
suite of models but the differences with the above result remain well within the uncertainty. 
All three models give equally good fits to the shape of the M band spectrum. An example 
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of such a fit is shown in Fig. 9 where the non-equilibrium spectrum corresponding to model 
C with log K zz = 4.3 is in much better agreement with the shape of the M band spectrum 
than the same model in equilibrium. The latter fails to produce the 4.665 /zm bump seen in 
the data that indicates the presence of CO. This corroborates our conclusion based on the 
fit of the 5-15 /xm spectrum that vertical mixing in the atmosphere of 2M ASS 0937 causes 
the chemistry to depart from equilibrium. 

By fitting absolute model fluxes to the data we obtain lower values of \ogK zz (Table 
2). If we rescale the observed fluxes by factors of 0.8-0.75 (less than the nominal calibration 
uncertainty), then models B through D can all be made to give logK zz (cm 2 /s) = 4.3 ± 0.3, 
and a consistent solution is found between fitting the shape of the M band spectrum and its 
flux level. This value of the eddy diffusion coefficient corresponds to a mixing time scale of 
~ 5 to 25 hours and a CO mole fraction of logX co = —5.10 ± 0.13. 



4-2.3. Photometry 



We attempted to fit the MKO V and M' magnitu des (iGolimowski et al.l 120041 ) and the 
Spitzer IRAC [3.6] and [4.5] fluxes (IPatten et al.l 120061 ) as in §4.1.2. In this case however, it 
was not possible to closely match the shape of the 3.5-5 /im peak with any of our models, 
not even with discarded model A (Table 2), for any value of K zz . In all cases, the modelled 
3.3/xni band of CH 4 is too deep, and the resultant V and [3.6] fluxes are too low (Fig. 10). 
Values of K zz that fit the [4.5] or the M' fluxes can be found but are different from each 
other. Using the model C parameters, we find \ogK zz (cm 2 /s) = 3.3 ± 2.2, but both the 
fit and x 2 are very poor, with the [3.6] flux falling several standard deviations below the 
IRAC measurement. The 3.3 /im band of CH 4 can be weakened by further decreasing the 
metallicity, or increasing K zz , but the latter to a value that is incompatible with the shape of 
the M band spectrum. Both of these possibilities appear unlikely, given the good agreement 
that we find with the spectroscopic data. A more complete line list for the CH 4 band would 
increase the opacity and thus go in the wrong direction. The discrepancy between our best 
fitting model and the [3.6] and V photometry might be attributable to the secondary effect of 
increasing the CH4 opacity and redistributing flux in the low opacity windows at wavelengths 
less than 1.7 /im into the 4-5 /im. 
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4-2.4- Summary: Optimal model 

In view of the difficulty in fitting the 3-5 fim photometry and the likelihood that the 
model spectra are impacted by an incomplete CH4 linelist, we again rely on the shape 
of the M band spectrum to determine the most likely value of K zz , and thus conclude 
that log-KT^ (cm 2 /s) = 4.3 ± 0.3. With the data in hand, we are unable to distinguish 
between models B, C, and D, as they give essentially equally good - or in the case of the 
photometry, equally poor - fits to the ensemble of data. Model C (T eff = 950 K, log g = 5.35, 
[M/H]=— 0.3) with logK zz = 4.3 as an example, provides an excellent fit to the red and 
near infrared spectra and to the Spitzer IRS spectrum as well (Fig. 11). It also gives a good 
fit to the shape of the M band spectrum (Fig. 9, although it underestimates the M' flux by 
1.4(7 (Fig. 10). This model also underestimates the [3.6], [4.5], and L' fluxes substantially, 
however. 

The chemical composition of the atmosphere as a function of depth is shown in Fig. 12. 
In the region of formation of the 4.7/xm band of CO (logT = 2.9 — 2.98), vertical trans- 
port increases the CO abundance by about two orders of magnitude. The NH3 abundance 
decreases by ~ 0.8 dex which accounts for the large effect seen in Fig. 8. The H 2 and 
CH 4 mole fractions are barely affected by non-equilibrium chemistry. As was found for 
Gl 570D, the nitrogen chemistry is quenched deep in the convection zone of the atmosphere 
(at logT = 3.38) and thus is unaffected by the value of K zz . On the other hand, the carbon 
chemistry is quenched much higher in the atmosphere, at logT = 3.12. 

5. Conclusions 

Low resolution M band spectra of Gl 570D and 2MASS 0937 clearly reveal the presence 
of photospheric CO via the CO absorption minimum between the P and R branches of the 

v — 1 band. These are the first spectroscopic detections of CO in late T dwarfs since the 

original discovery in the T7p dwarf Gliese 229B. The observed spectra are well-reproduced 
by model spectra with CO mole fractions of logX co = —4.10 ± 0.18 and —5.10 ± 0.13, 
respectively. 

Atmospheric models and synthetic spectra based on equilibrium chemistry predict that 
the CO abundance should be too low to form a detectable CO band and that instead the 
M band spectrum should be dominated by H2O features. Thus, the detection of prominent 
CO absorption in the spectra of these two T dwarfs is strong evidence for the presence of 
vertical mixing in the radiative part of their atmospheres. This interpretation is reinforced 
by independent evidence for depletions of NH 3 by a factor of 8 to 10 compared to the expec- 



-16- 



tations from chemical equilibrium calculations in both Gl 570D (jSaumon et all 120061 ) and 
2MASS 0937 (§4.2.1). The concordance of deviations of the NH 3 /N 2 and the CO/CH 4 ratios 
from chemical equilibrium in both objects can be explained consistently by a model including 
vertical mixing in the atmosphere and the kinetics of carbon and nitrogen chemistry. The 
present work firmly establishes this process in cool T dwarfs. 

The model spectra that include vertical mixing reproduce the entire observed SEDs 
of both dwarfs very well, with the exception of the photometry of the 4-5 /zm peak in 
2MASS 0937. Similar difficulties have been encountered with 2M ASS J04151954-0935066 
(T8) and 2MASS J12171110— 0311131 (T7.5) JSaumon et al-boOTl ). Given the overall success 
of our cloudless models at reproducing the SEDs of late T dwarfs, we think that this problem 
is related to well-known remaining flaws in the models. Since the 4-5 /im peak occurs in a 
low opacity window, it should be sensitive to the flux redistribution that we expect will occur 
in models that will benefit in the future from an improved line list for CH4 at wavelengths 
below 1.7 /im (see Fig. 7 for example). Our non-equilibrium spectra are currently computed 
as post-processing on a fixed atmospheric structure computed with abundances from chemical 
equilibrium. A self-consistent calculation of non-eq uilibrium models would als o redistribute 
the flux, but the effect should be relatively modest (IHubeny fc Burrows! 120071 ). 



The values of K zz and of r mix derived for Gl 570D and 2MASS 0937 depend directly 
on the choice of reaction pathway for the kinetics of the carbon chemistry and on the 
choice of qu enching sch e me. As in our previous work, we h ave adopted the "fast" reac- 
tion rate of lYung et all ( 19881 ) and the quenching scheme of ISmithl (119981 ) . The latter is 
physic ally and computationally more realistic than the s imple scheme that is generally ap- 
plied (IHubeny fc Burrows! 120071 : iFegley fc Lodderd I1994J ) . Its applicatio n to modeling the 
excess CO in Jupiter's troposphere t hen requires the fas t reaction rate of lYung et all (119881 ) 
to reproduce the spectroscopic data (IBezard et all 120021 ). We have used those same choices 
for the kinetics of carbon chemistry and the quen ching scheme in the past, so the present 
determinations of K zz and our p revious estimates ( Saumon et all 120071 ) are consistent. Us- 
ing the "slow" kinetic scheme of iPrinn &: Barshayl (119771 ) would result in significantly lower 
values of K zz and longer mixing time scales. 

The most direct way to determine the vertical mixing time scale r mix in the atmosphere 
with the data currently available is from the CO abundance derived by fitting models to 
the shape of the M band spectrum. The time scale is parametrized by the eddy diffusion 
coefficient, K zz = H 2 /r mix , for which we find \ogK zz (cm 2 /s) = 6.2 ± 0.7 for Gl 570D, 
which corresponds to a mixing time scale of 0.1 to 3 hours. While the uncertainty in the 
mixing time scale remains large, it is crudely comparable to the measured rotation period 
of P = 3.6 sin % hours, where i is the inclination of the axis of rotation to the line of sight 
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f lZapatero Osorio et al.lhood l For 2MASS 0937, \ogK zz (cm 2 /s) = 4.3 ± 0.3, implying a 



much longer mixing time scale of r mix = 5 t o 25 hours, which is also co mparable to the 



2-10 hour rotation periods of T dwarfs (jZapatero Osorio et all 120061 ). This suggests 



that the vertical mixing process is coupled to the rotation. The eddy diffusion coefficients 
of these two dwarfs differ by nearly two orders of magnitude, which seems a rather large 
variation for two T dwarfs whose physical properties are fairly similar: Alogg < 0.25, 
and AlogL = 0.22. Other properties in which there are significant differences between 
the two dwarfs in this study are the effective temperature, which differs by 100-150 K and 
metallicity, where A[M/H]=0.3 to 0.4. In the absence of a physical model for vertical mixing 
in brown dwarf atmospheres, it is difficult to appreciate the significance of these differences 
between Gl 570D and 2MASS 0937. In this context, a measure of the rotational velocity of 
2MASS 0937 would be of interest. 

Eddy diffusion coefficients (or bounds on their values) have been derived for the atmo- 
spheres of all of the solar system giants, often by comparing models with measurements of 
the vertical profiles of various photoche mical products. Other methods have been used as 



well. As reviewed by Moses et all (120041 ) . estimates for Jupiter's upper stratosphere cluster 
around \ogK zz (cm 2 /s)=6, comparable to or larger than the values we derive here for T 
dwarfs, despite Jupiter's four orders of magnitude smaller heat flow and nearly two orders 
of magnitude lower gravity. Constraints have even been placed o n eddy mixing in th e atmo- 



sphere of Neptune, where similar values of \ogK zz are inferred (IBishop et al.lll995l ). Such 
mixing in the radiatively-stable solar system stratospheres is often attributed to the propaga- 
tion and/or breaking of atmospheric waves. Bishop et al. review the literature in some detail 
and point out that K zz likely varies vertically through an atmosphere, increasing with falling 
density. Such a shortening of the mixing time scale with height would leave the quenching 
level and the resulting chemistry unchanged (Figs. 5 and 12). As such, the K zz values we 
derived from the observations are local values at the quenching level of the atmosphere, 
which is at both higher pressures and higher densities than the regions typically probed in 
giant planet atmospheres. The relatively similar mixing coefficients in such a great diversity 
of atmospheres is striking, even though the derived values apply to different atmospheric 
regions. Further investigation of the underlying mechanisms responsible for mixing in giant 
planet and brown dwarf atmospheres would be of interest and may also find applicability to 
the study of extrasolar giant planet stratospheres. 

Finally, we have derived the first set of well-constrained physical parameters for 2MASS 
0937 and have firmly established that this T6p dwarf has a subsolar metallicity and a mod- 
erately high gravity, as has long been suggested. Our analysis shows that, like other late 
T dwarfs, its mid-infrared and M band spectra can only be fit with models in which both 
the nitrogen and the carbon chemistry depart substantially from chemical equilibrium. This 
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reinforces the conjecture that vertical mixing and the associated departures from equilibrium 
chemistry are common, if not ubiquitous, among late T dwarfs. 

The data presented here were obtained at the Gemini Observatory, which is operated 
by the Association of Universities for Research in Astronomy, Inc., under a cooperative 
agreement with the NSF on behalf of the Gemini partnership: the National Science Founda- 
tion (United States), the Science and Technology Facilities Council (United Kingdom), the 
National Research Council (Canada), CONICYT (Chile), the Australian Research Council 
(Australia), Ministerio da Ciencia e Tecnologia (Brazil) and SECYT (Argentina). This work 
is also based in part on observations made with the Spitzer Space Telescope, which is oper- 
ated by the Jet Propulsion Laboratory, California Institute of Technology, under contract 
with NASA. Support for this work, part of the Spitzer Space Telescope Theoretical Research 
Program, was provided by NASA. We thank M.C Cushing for providing the Spitzer IRS 
spectra and A.J. Burgasser for the red spectra. This research has benefited from the M, 
L, and T dwarf compendium housed at DwarfArchives.org and maintained by Chris Gelino, 
Davy Kirkpatrick, and Adam Burgasser. 
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Fig. 1. — Spectra of T dwarfs 2MASS 0937 (upper panel) and Gl 570D (lower panel) near 
the CO 1-0 band center. Typical uncertainties are shown. Wavelengths of the 1-0 band lines 
and their relative strengths in optically thin gas at 750 K are shown in the upper panel. 
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Fig. 2. — Goodness of fit to the M band spectrum as a function of the eddy diffusion 
coefficient K zz for the nominal parameters of Gl 570D (T eff = 820 K, log ^ = 5.23, [M/H]=0). 
X 2 is defined in Eq. 1 (see text). The dotted curve is obtained when the model spectra 
are scaled to minimize x 2 j i- e - only the shape of the observed spectrum is fitted. The 
discontinuity in \ 2 near \o<gK zz (cm 2 /s) =7.2 occurs when the quenching of the carbon 
chemistry goes through the radiative/convective boundary in the atmosphere model. The 
solid curves are obtained when three absolute fluxes, corresponding to the nominal calibration 
(heavy line) and two that differ by ±30%, are included. 
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Fig. 3. — Fits to the M band spectrum of Gl 570D. The data are the black histogram and 
the typical pixel noise level (±lcr) is at lower right. The curves are model spectra plotted 
at i?=500 with T cS = 820 K, logg = 5.23, and [M/H]=0 that have been fit to the data with 
scaling factors adjusted to minimize the residuals. The red curve is the equilibrium model 
(K zz = 0), the blue curve is the non-equilibrium model shown in Fig. 4 (\ogK zz cm 2 /s = 4.5) 
and the green curve is the best-fitting model (logK zz cm 2 /s = 6.2). [See the electronic edition 
of the Journal for a color version of this figure.] 
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Fig. 4. — 3-6 spectral region of Gl 570D. Data are in black. Photometric measurements 



in th e IRAC [3.6], MKO L', and IRAC [4.5] bands (IGolimowski et all 12004! ; IPatten et al. 



20061 ) are shown by boxes whose widths approximate the bandpasses and heights indicate 
the ±1<7 uncertainty. Two models are shown at R = 200, both with T c g = 820K, log g = 5.23 
and [M/H]=0. The red curve is the equilibrium model (K zz = 0); the blue curve shows the 
model that best fits the photometry (logK zz (cm 2 /s)=4.5). Synthetic models fluxes in the 
filter bandpasses are shown by thick blue lines. The observed spectrum has been rescaled 
by a factor of 0.83 to match the model spectrum. [See the electronic edition of the Journal 
for a colorversion of this figure] 
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Fig. 5. — Chemical profile of the optimal model atmosphere for Gl 570D: T c g = 820 K, 
\ogg = 5.23, [M/H]=0 and \ogK zz (cm 2 /s)=6.2. Mole fractions are shown in equilibrium 
(K zz = 0, dashed curves) and out of equilibrium (\ogK zz (cm 2 /s)=6.2, solid curves) as 
a function of the temperature in the atmosphere (the top of the atmosphere is at left). 
Heavy black lines show the mixing time scale (r mix ) and the time scale for the destruction 
of CO ( t co , solid) and N2 (t N2 , dotted). The mixing time scale is nearly discontinuous 
where the atmosphere becomes convective (logT > 3.2). The 10-11 fim NH3 band is formed 
atlogT = 2.75 - -2.90 and the 4.7/xm band of CO is formed at logT = 2.8 - -2.9. [See the 
electronic edition of the Journal for a color version of this figure.} 
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Fig. 6. — T e fj and gravity for 2MASS 0937. Heavy black lines labeled with the mass in M Q 
are cloudless brown dwarf evolution tracks for [M/H]= —0.3. Isochrones (blue dotted lines) 
are labeled in Gyr. The nearly vertical lines (solid and dashed, red) are the locus of (T e s,g) 
points with \ogL/L Q = — 5.308±0.027. The solid dots show the four models of Table 2. The 
range of parameters for Gl 570D ([M/H]=0) is shown by the small box centered at 800 K 
( Saumon et al. 2006I ). [See the electronic edition of the Journal for a color version of this 
figure.] 
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Fig. 7. — Comparison of near infrared synthetic spectra corresponding to models A (red) 
and D (blue) in Table 2 to observed spectrum of 2MASS 0937 (thick black). Data are from 



Burgasser et al.l (]2006bl ). The models have a [M/H]= —0.3, do not include vertical mixing, 
and are plotted at R = 500. The main molecular absorbers are indicated. [See the electronic 
edition of the Journal for a color version of this figure] 
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Fig. 8. — Comparison of Spitzer IRS spectrum of 2MASS 0937 lCushing et al.l (120061 ). black 
histogram), best fitting model with equilibrium chemistry (model A, Table 2; red), and best 
fitting non-equilibrium model (model C, Table 2 with K zz = 10 4 cm 2 /s; in blue). The lower 
histogram is the noise spectrum. Flux calibration of the data is uncertain to ±3.7%. Model 
spectra are not renormalized to the data. The models have [M/H]= —0.3 and are plotted 
at R = 120. The main molecular absorbers are indicated. [See the electronic edition of the 
Journal for a color version of this figure.] 
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Fig. 9. — Fits to the M band spectrum of 2MASS 0937. Data are the black histogram; 
typical pixel noise level (±lcr) is shown at lower left. Models with T c r = 950 K, log g = 5.35, 
and [M/H]= —0.3, scaled to minimize residuals, are plotted at R = 500. The red curve is 
the equilibrium model (K zz = 0); the blue curve is the non-equilibrium model of Fig. 10 that 
best fits the shape of the spectrum (\ogK zz (cm 2 /s)=4.3). [^ee the electronic edition of the 
Journal for a color version of this figure.} 
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Fig. 10. — 3-6 /im spectral region of 2MASS 0937. Data are in b lack. Photometric mea 



surements in the I RAC [3.6], MKO U , IRAC [4.5], and MKO M' flGolimowski et alJ 12004 



Patten et al.l 120061 ) bands are shown by boxes whose widths approximate the bandpasses 
and heights indicate the ±1<t uncertainty. Two models are shown at R = 200, both 
with T eff = 950K, \ogg = 5.35 and [M/H]=— 0.3. The red curve shows the equilibrium 
model (K zz = 0) and the blue curve shows the model that best fits the M band spectrum 
(logK zz (cm 2 /s)=4.3). The synthetic models fluxes over the photom etric bandpasses are 
shown by thick blue lines. The distance is from IPerryman et al.l ( 19971 ). The observed spec- 
trum is scaled by a factor of 0.76 to match the non-equilibrium model spectrum. See text. 
[See the electronic edition of the Journal for a color version of this figure.] 



-32 - 



6 - 



£ 4 



* 2 - 








1.5 



2.5 



6 



£ 4 



* 2 - 








6 



8 10 
A (yum) 



12 



14 



Fig. 11. — Comparison of 2MASS 0937 spectrum (black histograms) with non-equilibrium 
model C (see Table 2) with T cS = 950 K, logg = 5.35, [M/H]=-0.3, with \ogK zz (cm 2 /s)=4.3 
(red curves). Spitzer IRS noise spectrum is the lower histogram in the bottom panel. Model 
spectra are plotted at a resolving power of R = 500 (upper panel) and R = 120 (lower panel) 
and are not rescaled to the data. Main molecular absorbers are indicated. [See the electronic 
edition of the Journal for a color version of this figure.} 
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Fig. 12. — Chemical profile of an optimal model atmosphere for 2MASS 0937 (model C of 
Table 2): T eff = 950 K, logg = 5.35, [M/H]= -0.3 and \ogK zz (cm 2 /s)=4.3. Mole fractions 
of H 2 0, CH 4 , N 2 , NH 3 and CO are shown in equilibrium [K zz = 0, dashed curves) and 
out of equilibrium (logi^ (cm 2 /s)=4.3, solid curves) as a function of temperature. Heavy 
black lines show the mixing time scale (T mix ) and the time scale for the destruction of CO 
(t co , solid) and N 2 (r Na , dotted). The mixing time scale is nearly discontinuous where 
the atmosphere becomes convective (logT > 3.25). The 10-11 /im NH3 band is formed at 

logT = 2.80 2.95 and the 4.7 /xm band of CO is formed at logT = 2.90 2.98. [See 

the electronic edition of the Journal for a color version of this figure.} 
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Table 1. Observing Log 



UT Date 


Source Name 


Exposure (min) 


Calib. star 


20040429 


Gl 570D 


64 


HIP 78585 


20040430 


Gl 570D 


64 


HIP 78585 


20041230 


2MASS 0937 


90 


HIP 52959 


20050116 


2MASS 0937 


84 


HIP 52959 


20050126 


2MASS 0937 


108 


HIP 52959 


20050215 


2MASS 0937 


84 


HIP 40843 


20050222 


Gl 570D 


105 


HIP 67945 


20050301 


Gl 570D 


96 


HIP 77610 


20050302 


Gl 570D 


63 


HIP 67945 



Table 2. Range of Physical Parameters for 2MASS J09373487+2931409 a b 



Model 


T e ff 


logs' 


logL/L 


Mass 


Radius 


Age 


\ogK zz c 




(K) 


(cm/s 2 ) 




(Mj) 


(Re) 


(Gyr) 


(cm 2 /s) 


A 


865 


4.86 


-5.296 


27 


0.0984 


1.0 




B 


923 


5.20 


-5.304 


45 


0.0870 


3.0 


3.2 ±0.10 


C 


950 


5.35 


-5.308 


57 


0.0818 


5.0 


3.0 ±0.11 


D 


974 


5.47 


-5.312 


69 


0.0778 


10.0 


2.7 ±0.14 



a [M/H]=-0.3 

b The la uncertainty in the luminosity is AlogL/L = 0.027, corresponding 
to AT cff = 15 K and AR/R Q = 0.001 at constant gravity. 



c For the nominal flux calibration of the data. 



